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THEORETICAL ELASTIC STRESS DISTRIBUTIONS ARISING FROM DISCON- 
TINUITIES AND EDGE LOADS IN SEVERAL SHELL-TYPE STRUCTURES 

By Robert H. Johns and Thomas W. Orange 


SUMMARY 

The deformation and complete stress distribution 
are determined for each of the following edge-loaded 
thin shells of revolution: (1) a right circular cylin- 
der , ( 2 ) a frustum of a right circular cone , and 
(3) a portion of a sphere . The locations of the 
maximum circumferential and meridional stresses 
on both the inner and outer surfaces are also found. 
The basic equations for the above were selected from 
the published literature on the subject and expanded 
to produce resultant-stress equations in closed form 
where practicable to do so. Equations are also devel- 
oped for the discontinuity shear force and bending 
moment at each of the following junctions: ( 1 ) axial 
change of thickness in a circular cylinder , (2) axial 
change of thickness in a cone , (S) change of thickness 
in a portion of a sphere , (4) a cylinder and a cone, 
(5) a cylinder and a portion of a sphere, (6) a 
cylinder and a flat head , and (7) a cone and a 
portion of a sphere . 

INTRODUCTION 

Weight considerations in space-flight structures 
require loading or stressing of the structures very 
nearly to their maximum capabilities. It is 
therefore important that the operating stresses be 
known with a high degree of accuracy. With 
this knowledge it is possible to obtain the best 
ratio of structural weight to gross weight and to 
ensure the structural integrity of the vehicle. 

Shell structures offer excellent weight and fabri- 
cation characteristics for use in missile and space 
structures. It is the purpose of this report to 
present methods of analysis for several problems 
encountered in this type of structure. The mem- 
brane stresses produced by the pressurization of 


such shells are usually easily computed. How- 
ever, the forces and stresses involved in the dis- 
continuity regions are not so readily determined, 
and the published techniques of solution are in 
many cases of such a nature as to preclude then' 
use by the design engineer. There are four basic 
difficulties involved. First, for some problems no 
solution exists in the published literature. Sec- 
ond, many of the published solutions involve 
mathematical complexities beyond the background 
of the average designer. Third, a large propor- 
tion of the solutions, as presented, are not carried 
to the point where they can be used directly for 
determining stress distributions. Fourth, many 
of the solutions, depending on the assumptions 
involved, require calculations of an extensive and 
tedious nature. 

In general, the more rigorous solutions are based 
on few assumptions and lead to difficult analyses 
which are of no practical use to the designer. It 
is necessary to sacrifice some degree of accuracy to 
obtain a solution in a reasonably direct mathe- 
matical manner. The degree of accuracy that is 
required and the degree of complexity that is con- 
sidered acceptable are factors which must be 
weighed in selecting a method of analysis lor any 
shell problem. 

An attempt was therefore made in the prepara- 
tion of this paper to survey the literature to sort 
out the more practical typos of solutions. It was 
necessary to expand some of the solutions exten- 
sively to obtain expresssions for all the internal 
forces. From these, complete stress distributions 
and the maximum stresses and their locations 
were determined. The selected solutions, as 
modified and extended, are compiled in this report 
in a form which can be readily used by the average 

l 
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engineer. These solutions include treatment of 
the following problems : 

(A) Deflection, rotation, and internal stresses 
in the following edge-loaded shells of 
revolution : 

(1) Right circular cylinder 

(2) Frustum of a right circular cone 

(3) Portion of a sphere 

(B) Discontinuity shear and moment at the 
following shell junctions : 

(1) Axial change of thickness in a circular 
cylinder 

(2) Axial change of thickness in a cone 

(3) Change of thickness in a sphere 

(4) Junction of a cylinder and a cone 

(5) Junction of a cylinder and a portion of a 
sphere 

(6) Junction of a cylinder and a flat head 

(7) Junction of a cone and a portion of a sphere 

The analyses for the changes in thickness are 

intended as an investigation of the edge effect of 
chemically milled or machined surfaces where an 
abrupt change of thickness symmetrical about the 
middle surface occurs. Several special cases can 
be deduced from those just given. The equations 
for the edge-loaded portion of a sphere can be 
specialized to give the case of an edge-loaded 
hemisphere. The solution for the junction of a 
cylinder and a portion of a sphere will yield the 
special case of a cylinder with a hemispherical 
dome. The solution of a cylinder with fixed ends 
can be obtained from the case of the cylinder with 
a flat head. The solution of a junction of a cone 
and a portion of a sphere has the special case in 
which the cone and spherical shell have a common 
tangent to their meridians at the junction. This 
solution combined with the solution for a cylinder 
with a hemispherical dome will yield the solution 
of a toriconical head in which the torus is a portion 
of a sphere. 

Only the final equations and parameters neces- 
sary for the solution of any of the cases are given 
in the body of the report. The derivations, to- 
gether with intermediate equations and some dis- 
cussion, are presented as appendixes A to L. 

SYMBOLS 

a radius of cylinder, in. 

c thickness ratio 

D Eh? j\2(\—v 2 ) , lb-in. 

E modulus of elasticity, psi 

77 radial shear force in wall of shell acting on a 


plane perpendicular to axis of revolution, 
lb/in. 

h thickness of shell wall, in. 

M bending moment in wall of shell, in.-lb/in. 

m 12(1 — v 2 ) 

N uniform normal force in wall of shell, Ib/in- 
p internal pressure, lb/sq in. 

Q shear force perpendicular to wall of shell, 
lb/in, 

U radius of spherical shell, in. 
r radius of parallel circle, in. 

V angle of rotation of a tangent to a meridian, 
radians 

w deflection perpendicular to axis of revolu- 

tion, in. 

x distance along meridian from edge of cylin- 
der, in. 

y distance along meridian from apex of cone, 
in. 

a half-angle of cone or portion of sphere, 
radians 

ft </3(l-v 2 )/a 2 h 2 , in. -1 

& edge-deflection influence coefficient 

0 ( ) e~ ( ) cos ( ) 

X ^3(1 -v 2 )R 2 /h 2 

fx 12(1 — v 2 )lh 2 tan 2 a, in.~ I/2 

v Poisson’s ratio 

£ 2m Vy 

c t normal stress, psi 

r shear stress, psi 

$ ( ) €“ ( } [cos( ) + sin ( )] 

<p angular measure in plane containing merid- 

ian and axis of revolution, radians 
^ ( ) e~ ( } [cos( )— sin ( )] 
ft ( ) e~ () sin( ) 

o) edge-rotation influence coefficient 

Subscripts: 
c cylinder 

f flat head 

i , j dummy indices 

k cone 

m meridional 

s spherical shell 

x meridional direction on cylinder 

y meridional direction on cone 

6 circumferential 

r shear 

<p meridional direction on spherical shell 

0 junction 

1,2 different thicknesses at change of thickness 
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Because some subscripts are used only once or 
twice, they are, for purposes of brevity, defined 
where they occur and not listed here. This is 
particularly true for many parameters appearing 
in the cone analyses. Also, the usage of a number 
of subscripts is considered self-evident, and no 
amplification of their meaning is given. 

METHOD OF ANALYSIS 

In the design of pressure vessels, certain regions 
frequently exist where continuity of the structure 
cannot be satisfied by membrane forces alone. 
Such regions are known as “ discontinuity” areas. 
The discontinuity forces which are induced to 
make the rotations and deflections of the walls 
continuous are usually of a local nature, but they 
may considerably alter the stress distributions in 
the regions where they occur. It is assumed in 
this report that, where more than one discontinuity 
is present, the distance between them is sufficient 
that each does not noticeably affect the discon- 
tinuity shear and moment of the other. 

All shells are rotationally symmetric and loaded 
by internal pressure. The stresses due to the 
weight of the structure are not considered. These 
are usually much smaller than those due to inter- 
nal pressure and may frequently be neglected 
without introducing noticeable error. The 
stresses due to dead weight, supports, concen- 
trated loads, or other such conditions can be 
superimposed upon those presented here when such 
conditions exist and are significant. 

The usual method of determining the discon- 
tinuity forces is to imagine the shell to be physi- 
cally separated at the discontinuity. The edges 
of the two components will, in general, rotate and 
deflect different amounts if the membrane stresses 
alone are considered. Deflections and rotations of 
the edges of the components can be found from 
conventional membrane analyses. To preserve 
continuity of rotation and deflection in the actual 
structure, a discontinuity shear and moment must 
be present on the edge of each component. Equi- 
librium of forces at the cross section requires that 
the shear and moment on the edge of one compo- 
nent be equal and opp site to those on the .mating 
edge of the other component. Having expressions 
for the edge rotation and deflection of each 
component due to edge shear and moment, as 
well as to internal pressure, it is possible to 
write two equations expressing the equality 


of deflection and rotation at the discontinuity. 
These two equations can be solved for the un- 
known discontinuity shear and moment. Once 
these have been determined, it. is possible to 
find the rotation, deflection, and internal forces 
for any element of the shell. 

The edge-loaded shells and the regions of the 
shell on either side of the discontinuity are as- 
sumed to have constant thicknesses. The sur- 
face described by the revolution of a meridian 
midway through the thickness is called the middle 
surface. The middle surface of the shell is as- 
sumed to be continuous from one region to another 
across the discontinuity. Hence, no additional 
moments are induced by mismatching of the 
effective lines of action of the meridional forces in 
the two regions. A discussion and analysis for 
including the effect of nonconcurrence of the 
middle surfaces at the junction is presented in 
appendix L, but it is not used in the bod3 T of the 
report. The shear forces acting on sections made 
by planes containing the axis of revolution and 
intersecting the shell are zero because of axial 
symmetry. Thus the circumferential or hoop 
stresses acting on these meridional planes (planes 
containing the axis of revolution and a meridian) 
are principal stresses. 

The normal stresses perpendicular to the middle 
surface of the shell are neglected, since they are 
usually much smaller than those in the circumfer- 
ential and meridional directions. These radial 
stresses which are neglected vary in magnitude 
from the value of the pressure on the inside 
surface of the shell to zero on the outside. Con- 
sequently, since the radial stress is neglected, 
a biaxial stress state is assumed. The direct 
stresses are assumed uniform throughout the 
thickness. The bending stresses are assumed to 
vary linearly through the thickness from zero at 
the middle surface to maximum values at the 
inner and outer surfaces. 

The shearing stress varies parabolically through 
the thickness with the maximum occurring at the 
middle' surface and decreases to zero at the inner 
and outer surfaces. As mentioned before, no 
shear stresses act on meridional planes. The 
meridional stresses, that is, stresses acting on 
planes normal to the meridian, are very nearly 
equal to principal stresses. This is because the 
shear stresses are usually much smaller than the 
normal stresses. In fact, the meridional stresses 
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on the inner and outer surfaces are principal 
stresses because the shear stresses are zero at these 
surfaces. 

Because of the nature of the bending stresses, 
the maximum and minimum stresses will be found 
on the inner and outer surfaces. When solving 
for the locations of the maximum stresses, and 
from these the stresses themselves, the usual 
theory of maximums and minimums is employed. 
Since the equations for stress are exponentially 
decaying sine functions, more than one solution 
will be found for the possible location. However, 
only the smallest positive value of the independent 
variable so obtained is significant. 

The maximum stress sometimes occurs at the 
loaded edge of the shell when the slope of the com- 
bined stress curve has a nonzero value. There- 
fore, when solving for maximum stresses, it is 
always necessary to check the loaded edge of the 
shell in addition to the location of zero slope 
nearest to the loaded edge. Because of the bend- 
ing stresses, both the inner and outer surfaces must 
be checked for the maximum stress. By checking 
both surfaces, both the maximum and minimum 
stresses will be obtained. 

The material is assumed to be homogeneous 
and isotropic and to obey Hooke’s law. The 
results are applicable only to stresses within the 
elastic region for thin-walled pressure vessels. 
Also, small deflections are assumed throughout 
the report in the derivation of the basic equations. 

SUMMARY OF SOLUTIONS 

The summaries are arranged here in the order 
Tn which they would ordinarily be used. First, 
the discontinuity shear and moment are deter- 
mined for the particular shell junction being con- 
sidered. Having determined these, it is only 
necessary to substitute them into the corre- 
sponding edge-loaded shell equations to deter- 
mine the stress distributions. The arrangement 
of the summaries is not the same as the order in 
which they are derived in the appendixes. This 
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is because some of the equations derived for the 
edge-loaded shells were necessary for the solu- 
tion of the junction problems. Consequently, 
the edge-loaded shell solutions appear first in the 
appendixes. 

DETERMINATION OF DISCONTINUITY SHEAR 
FORCES AND BENDING MOMENTS 


Shear and moment at an axial change of thick- 
ness in a circular cylinder. — Referring to figure 1 
for sign convention, the discontinuity shear and 
moment for this case are 


and 


where 


Qo- 


.l/o= 


2-v 

r (c— i)(c 5/:! +]) 

' 20, 

L(c J +l)*-f2c 3 ^+l). 

2-v 

r (e-l)(e*-l) ' 

4/s; 

L(e*+l)*+2<"*(c+l). 


p (1) 


V ( 2 ) 



(3) 


and the subscript 1 refers to the cylinder with 
thickness h\. Because of the sign convention 
chosen in figure 1, the sign for Q 0 must be changed 
when solving for the stresses in the cylinder of 
thickness h 2 . 


t 


jhL 


a-, rM- rO n 
o Vs' 0 \/ u 






^7 


h 2 P 




\V pa '• 
2 


1 






Middle 
| surface 
o 


4 


Axis of 
revolution 


Figure 1. — Change of thickness in cylinder. 


Shear and moment at an axial change of thick- 
ness in a cone. — -With the sign convention as 
shown in figure 2, the discontinuity shear and 
moment for this case are 


and 


jj (< 5 i.a / 0 5 2 i ji/ 0 )(wi iP w 2 . p ) — (81 tP — $2, *>)(wi,Ar 0 — w 2 j a/ q ) 


($1.J5T 0 “ $2. JT 0 ) — «2.Jir 0 )“*($l,JI# 0 — 5 2 ,Af 0 )(wi./f 0 — P 


(4) 


(5) 
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The subscripts 1 and 2 refer to the regions of 
thicknesses h x and h 2 , respectively. The expres- 
sions for the edge influence coefficients are ob- 
tained from their more general counterparts pre- 
sented in the section “Frustum of a cone loaded 
by edge shear and moment and internal pressure .’ 1 
For the problem being considered here, yi=y 2 —y^ 
and r 1 =r 2 =r 0 , but ^ and £ are not the same for 
the two edges at the junction, because h x does not 
equal h 2 . Therefore, 


fi=2jxrvyo 

f 2 =2^2\yo 


m 


MU 


M2 = -7 



h >2 tan of 


( 6 ) 


The parameters X b £i, «i, and are associated 
with the region of thickness h x ; likewise, X 2 , £>, 
< t) 2f and 5 2 are associated with h 2 . 

The expressions for the edge influence coeffi- 
cient are 


. rn z j 2 r 0 
EhfVhi cosa fil 


_m 2 r 0 


'*» Eh\ 


«i, 


_mVotana 0 3(l+i»)tan*a / r o p 3r 0 tanj* 

p- 2 Eh\ i+ mE \2AiCos« 1 2Eh 1 cosa 


v\,M n ma **2 


Eh\ 

_ mr o l2r 0 cos a 

5l "° — EhS ~hT~ 


Si. 


in a j r Q 0 m 2 r 0 tan 2 « „ , 
Ehi V 2Ai cos a 0/1 2i * 


mro sin 




8(1 — r)Z? Eh\ cos a 


_ m 3 j 


2r f> 


A 2 cos a 




m 2 r 0 


mVotana^ 8(1 +r) tan 2 a / r 0 

4 " 


0)2,1 


2 Ehi 


: 0.- 




3r 0 tan a 
2£ r /t 2 cosa 


, _Z^o 

£7tr 25 


„ mr„ /2r 0 cos a „ 

"u 9 


, mr 2 sin a / /•„ 0 

Eh, V2A 4 cosa 6 


m 2 r (l tan 2 




S(l — p)/? ° COS a J 


( 7 ) 


50C835— r,l i 
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where 


n 1 = 


bG 


Qo - 


y*l{C+2vG) 

-A 


C+2vG 

£1 

\V(0+2vG) 


12 


£2 GjC 


12 


4 ^(C K +2vG K ) 

— A K 

vGk 

-^b k 


^ ( 8 ) 


4 +Q k 


12a 


V2 (C k +2vGk) . 


and 


A=h(ber'<£i behti— bei'^i ber 2 £ i) 

B— (ber'£i) 2 + (bei^i)* 

C =£ i (ber 2 £i ber^ + bei^i bei^i) 
G=(ber 2 h)*+(bei 2 hY 

l (9) 

A K =£ 2 (ker' 2 Z 2 kei£ 2 —kei r £ 2 kerfo) 

B K =(ker’&)*+(kei ’&) 2 

C K =Z 2 (ker 2 £ 2 ke r' 2 £ 2 +kei 2 ^ kei^ 2 ) 

G k = (ker 26) 2 + (ke ^ 

Figure 3, taken from reference 1, can also be used 
to determine graphically the 12 functions for cer- 
tain ranges of £. If the analytical expressions are 



used, reference 2 can be used to find values for the 
Bessel-Kelvin functions and their derivatives up 
to arguments of 107.5. For large arguments, the 
following asymptotic expressions may be used : 


12 1 


VS 1.855 

io£ e 


n 2 


VS 3.730 

10£ {* 


^3 


0.84853 1.865 

4 e 


Q 4 «l 

n,=i 

n#=i 


VS 1.855 ■) 

mos j 2 

VS 3.730 

^lOi £ 2 

0.84853 1.865 

+ i e ) 
( 10 ) 
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The desired discontinuity shear and moment are 
now found by use of equations (4) and (5), respec- 
tively. Essentially all of the preceding analysis 
is taken directly from reference 1. The nomen- 
clature is essentially the same as in reference 1 
except for some differences due to changes in sign 
convention. 


Shear and moment at a change of thickness in 
a portion of a sphere— The discontinuity shear 
and moment for this case are 


TJ (1 —v)R 

r a-c)(i+c«) i 

0 o\ 

2Xi sin a 

|_(l+c 2 ) 2 +2c 3/2 (l+c)J 

and 


w (i— ‘ 

- (1— c)(l— c 2 ) -] 

-VJ° 4X 2 [ 

_(l-f c 2 ) 2 +2c 3/2 (l+c)J 

where 



P 


P 



(id 


( 12 ) 


(13) 


used in figure 4, the sign for H 0 must be changed 
when solving for the stresses in the region of thick- 
ness A 2 . Also, the angle (ir— a) must be substi- 
tuted for a in the stress equations for thickness ho. 



Figure 4. — Change of thickness in sphere. 


and the sign convention is as shown in figure 4. Shear and moment at the junction of a cylinder 
The subscript 1 refers to the portion of the sphere and a cone. — In this case the discontinuity shear 

with thickness h\. Because of the sign convention and moment are 


and 


° $c,H 0 ) (w*. Mq—UcMq) ~~ °c,M 0 J \^k.H 0 ~ - 

_ (<$*, p ~~ &c. p) (g*, H 0 —Ue, Hq) ~ y (frt, H 0 ~ &c. H 0 ) 

° (Sk,H Q — &c.H Q )(Uk,M 0 — C0 e ,3/ 0 )~ (5* t A/ 0 — S CtAfo )(o3 k ,H Q — 


(14) 

(15) 



Figure 5. — Junction of cone and cylinder. 


where the subscripts c and k refer to the cylinder 
and cone, respectively. Most of the edge in- 
fluence coefficients appearing in the preceding 
equations are evaluated from the expressions 
found in the summaries for the edge-loaded 
cylinder and edge-loaded cone (eqs. (30) and 
(38)). In addition, from membrane theory 





The sign convention for II 0 and d/ 0 is given in 
figure 5. 




(a) Geometry and sign convention. 

(b) a <90°. * 

(c) a = 90°. 

(d) or>90°. 

Figure 6, — Junction of cylinder and portion of sphere. 
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Shear and moment at the junction of a cylinder 6(a) for the sign convention, the discontinuity 
and a portion of a sphere. — Referring to figure shear and moment for this shell junction are 


Z/ 0 = 


ffi — - — (2 — v)cl [1 +c 2 >c sinal — — [l+c 2 +2c 2 ^c sin a] cos a 
J \_sm a J A* 

2r 2 

(l+c 2 ) 2 + , . - (1+c sin a) 


\C sin or 


pR 


and 


M o = 


f >-] MHt f saa) a “ A e* 


(1+C 2 ) 2 


2c 2 

\C sin a 


(1+c sin a) 


(17) 


(18) 


where the subscript c refers to the cylinder and 
$ to the portion of a sphere and 



(19) 


Equations (17) and (18) are true whether a is 
greater than or less than 90° (see fig. 6). For 
the particular case a =90°, or R=a , the solution 
is obtained for a cylinder with a hemispherical 
dome. Note that, when solving for the stresses and 
displacements in the spherical portion of the shell, 


pR 

Hq+— COS a 
J* 


is substituted for H 0 in equations 


(42) to (52) (see fig. 6(a)) for this case. 

Shear and moment at the junction of a cylinder 
and a flat head. — The discontinuity shear and 
moment for this case are 


/ c 3 \c+2(2—p)c 3 X c +2(2 — p)(1+p) \a 

\ 2c 3 Xc+[(l— i')c 4 +(l+i')]Ac+(l““*' 2 )c J 4 

( 20 ) 


f 2c 3 X 3 + ( 1 — y)c 4 Xc +2 (2 — p) ( 1 + y) 

* / o “ \ 2c 3 X2 + [(l-.)c 4 +(l+^)]X c +(l-^)c J 8X C 1 

( 21 ) 


where 



( 22 ) 


and the subscripts c and j refer to the cylinder 
and flat head, respectively. The sign convention 
is given in figure 7. If the thickness of the 
flat plate is very large in comparison with the 
thickness of the cylinder (h f ^> h c ), c approaches 
zero and the solution for a cylinder with fixed 
ends free to expand axially is obtained. 



Figure 7. — Junction of cylinder and flat head. 
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(e) (f) (g) 


(a) Geometry and sign convention, 

(b) <*,>! —a t . (c) a,= (d) a.<|— a k . 

(e) o,>|+a t . (f) a, = |+ot. (g) «,<|+or t . 

Figure 8. — Junction of cone and portion of sphere. 
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Shear and moment at the junction of a cone 
and a portion of a sphere. — Using the sign con- 


vention shown in figure 8, the discontinuity 
shear and moment for this case are 


and 


M 0 =- 


{ tt 

[w«,J? 0 (5*,.V 0 — 5,,H 0 (uk.M 0 — "s..v 0 )l ~2 C0S 

-f- (5* , p — 8, t p ) ((ilk, M 0 (&k,i r 0 5j.i/ 0 ) (“*. p w *. ») jP 

S,.M 0 )(uk.no— — 




(5t.4f 0 — Kit 0 ) (o>k.H 0 —W.,H 0 )—(&k,H' ) —5,.H ( )(Wk,M 0 —<il,.M a ) 


(23) 


(24) 


where the subscripts k and s refer to the cone and 
the portion of a sphere, respectively. From 
membrane theory 


«».p=0 


and 


K 


l-p R] 
2 Eh, 


sin a. 


(25) 


The remaining edge influence coefficients are ob- 
tained from the summaries for the edge-loaded 
cone and portion of a sphere (eqs. (38) and (50)) 
with changes in the signs for and due to 
the change in sign convention for V. Also the 
signs for and 8>c.h, must be changed because 
of the change in sign convention from II, to H 0 . 
Because of the sign convention chosen in figure 8, 
the sign for II 0 must also be changed when it is 
substituted for H x in the equations for the edge- 
loaded cone to find the stresses (see also fig. 11). 

pR 

Note that Ho+^r cos a s is substituted for H 0 in 
& 

equations (42) to (52) when solving for the stresses 
and displacements (see fig. 6(a)). The cone and 
portion of a sphere need not have a common 
tangent to their meridians at the junction. For 
the special case where the cone and spherical shell 
are tangent and the spherical shell is tangent to a 
cylinder of the same radius, the solution is obtained 
for a tori conical head that has a torus which is a 
section of a sphere. The various possible sym- 
metrical junctions of conical frustums and 
spherical shells are shown in figures 8(b) to (g). 


DETERMINATION OF STRESSES AND DEFORMATIONS IN 
EDGE-LOADED SHELLS 

Right circular cylinder loaded by edge shear and 
moment. — The following equations are applicable 
at any distance x from the loaded edge of the 
cylinder. The type and location of stress, as well 
as the rotation and deflection, are given by equa- 
tions (26) to (29). Where ± or =F signs occur, the 
upper sign refers to the stresses on the inner 
surface and the lower sign to the stresses on the 
outer surface. See figure 9 for the sign convention 
and figure 10 for curves of the functions 0, 
and 12. Values of these functions are also tabulated 
in references 3 and 4. 

Meridional stress : 


* I =± ![M 0 <W+i<?oSW 


(26) 


Circumferential stress : 


=[ 20 * l *<fix) ± *(jfa)] Af,+[20 jjeote) 


6v 

W 


(27) 


Rotation of meridian: 

V= WD hW) + <?o*G&0] (23) 

Outward displacement : 

[™) + f? 00 ^)] (29) 
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+y 

w 



(a) Geometry and loading. 

(b) Internal forces and moments. 

Figure 9. — Cylinder. 


At the loaded edge, these equations reduce to 


£ 

N 

1 

o 

II 

H- 

6 M 0 

(26a) 

A 2 


)M O +20^Q O 

(27a) 


Fo — 2 fpD ( 2 ^ M o+Qo) =c * > M 0 Mo J r ( * ) Q 0 Q 0 


W °~2^D G^o+Go) ~^m 0 M 0 -\-5q q Q 0 


(28a) 


(29a) 


where 


wj# « = w> 


i 

“°o-o pD 

So,=7 1 


0 2 pD °° 2fi 3 D J 


Tlie locations of the peak stresses are given by 
the following equations : 

Meridional stress: 


Circumferential stress : 


(31) 


Xe— - arc tan 
P 


T rn 2 M 0 + (W Gy) 
l2^Ai 0 -f- (± w 2 +Gy)(? (l _ 


(32) 


Frustum of a cone loaded by edge shear and 
moment and internal pressure. — The following 
equations are applicable at any distance y from 
the apex of the conical frustum where yi'>y>y 2 
(see fig. 11). The type and location of stress, as 



(30) 


Figure 10. — Functions used in analysis of cylindrical 
and spherical shells. 





(b) 

(;i) Geometry and loading. (b) Internal forces and moments. 


FrouHK 11. ---Frustum of conr. 

well as the rotation and deflection, are given by of the Bessel-Kelvin functions as tabulated in 

equations (33) to (36). Where i or T signs reference 2. The primes denote differentiation 

occur, the upper and lower signs refer to the with respect to £. Solutions for the constants of 

stresses acting on the inner and outer surfaces, re- integration C 2 , ( and (\, as well as other 

spectively. See figure 11 for the sign convention parameters appearing in the following equations, 

and loading. The solutions arc given in terms arc given at the end of this summary. 

596835—91 — — 3 
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Meridional stress: 


a * = hy\ ^ nP + 

+^’2 jTiofi"; (£ 6w2^-t-2«'6e/j|)J 

+ ^ :* t T — (s^ v^ + 2i'/v'y.,|)J 

. tan a /« , 3 tan ct\ 

+— ( s ^ u ^)'' 


( 33 ) 


Circumferential stress : 


=^{ r ’> [I (26^+^*)] 

+ C 2 [^| bei’^±-^ (2ber 2 Z+vZber,Z)~^ 
+<~' 3 ker'^T — {2kcut+vtkei 2 £)~^ 

+ (\ ±^jjj (2&er 2 $+»»|&er£|)J^ 

, . /v , 3 tan a\ 

+tana ^ ±iT _^ 


(34) 


Rotation of meridian: 


rn 1 

T 7 =-gp (Ci bei 2 £—C 2 ber 2 Z+C 3 ke i 2 £ 


(35) 


Outward displacement : 

^ = j?/r {/'' (I htr 'it— vber *i) 

+ c* Q 2 bei'jt—pbe.L^+Ci ker' 2 £—vker,{^ 

+r. (‘ )]+g£s (,_|) 


V 

m) 


These equations reduce to the following equa- 
tions at the loaded edges, where the subscripts ] 
and 2 refer to the edges and ?/ 2 , respectively, 

Wj sin a cos a sin a TT . 6 w /IJO x 
ff », =- 2A 77 h~ 7/1 ± p M ‘ (3 - ia) 


iji sm a cos a sin a 


2h 


tj , 64 /, 
p h Hi± h* 


(33b) 


f\ 


V /~T 1 A /"r 


’> 4 1-r \ C+2W7 

/ 

7/1 tan a 


1 — — v * ^ t ? ^ 

±6^.4-^ B-vC 


sm- 4 a m L 


h 


Sill a 


1- 


±6L 5f-^-# £ -' r 


C+2rG 


C+2vG 
Hx 


P 


1 


:2g\ 


„ -j.m i A± v C±2G } 

+ P \ C+2r(? ) Ml 


(34a) 


3 tan 2 a /I w7 2 ^4a; -F vG k -^-^Gk . 

2 ~4 \~V \6 C k +2pG k ) P 


1/2 tan a f j ^sin 2 a 


±6 


1-V 


m 


4 if O „^nr ^ 

2 — 7) n K — ^ 


! 2 
1 — v 2 


C k -\-2vGk 


. / ±6 
sm a I m 


J If 77 N 

2 — 77 &K — v ' K 


1 


Ck-\-2vGk 

Fi=wi t ji/|Mi +a>i f Hi^\ 4"Wi f vV 
T 2 = 0)2, 4/ 2 4C + 0) 2 , H'fll ^2, p P 

W >2 = & 2 , M 2 M 2 + §2 , H 2 H 2 + & 2 , V P 


6 [ m2 AK±vC K ±2G K A 

H2+ h\ C k +2vG k ) Mi 


(34b) 

(35a) 

(35b) 

(36a) 

(36b) 
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where 


r,=- 


a=- 


r 3 =- 


r 4 =- 


(^i2\?/iA sin a+^ sin 2 a tan 

(£, ber’ 2 Zi+ 2v ber 2 ^)+2m 2 y l bei 2 $\ 

~ w ^A 2 tan 2 a“| 

8(l-r) J 

(ihyih sin a+| pAj/f sin 2 a tan aj 

h(C+2vG) 

(£i bei'fo+2v bei 2 h)—2m 2 yi ber-& 

", r ph 2 tan 2 af] 

_ 1 8(1-0 J 

(Hzyjl sin a+^/>/i ?/2 sin 2 a tan 

h(C+2vG) 

) (£2 ker’£ 2 +2v ker 2 Q +2 m% kei 2 h 

r, r M 2 tan 2 oT| 

L M ’ 8(1— f) J 

(H 2 y*h sin a+~ 7 >A.y| sin 2 a tan a ^ 

k(C K +2 v G K ) 

) (fc kei’£i+2v kei 2 £ 3 ) —2m 2 y- 2 ker.& 2 

r ph 2 tanV] 

l Mi 8(1-0 J 

h{C K +2vG K ) 


^ (37b) 


7c) 


and 


™l r o 

Eh 2 \hc os a 1 


— m ri o 
^ ”£A 2 


0)! 


m 2 rf tan a 0 t 3 (1 + p) tan 2 a j n ^ 3n tan a 
2Eh^~ mi? \ 2A cos a 1 2Ek cos a 


<5l Jfc 


_ m*ri 


* EH 1 ^1 cos a Q-i 

A’A\ A 


mr] sin a / ri ^ m 2 ri tan 2 ^ ^ 


...('-O' 1 


^ l p Eh \2h cos a 3 S(l — »0i? * Eh cos a 


m 3 / 2r 2 0 

( * >2 > M 2~Eh 2 \ h cos a 4 


m*r.. 




^ Eh 2 
m' 


r\ tan a Q __3 (l+») tan 2 a j 3r 2 tan cy 

2i?A 2 ' J mi? \2Acosa 4 2/?A cos a 


m z r> 


/t7( 2 ’ fl * 


mr 2 /2r, cos a > 

5i «>~Krv — r“ * 


. mr* sin a / 0 

5 -'^ £7t \2Acosa 8 


m 2 r 2 tan 2 « ^ . N 

^ S2 5+^j 


(-!)'• 


8(1-0® ®A cos a 


(3S) 
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The remaining parameters necessary for the solu- 
tion of the cone are given in equations (8) and (9). 
Equation (10) can also be used where appropriate. 
As previously mentioned, the cone analysis was 
taken directly from reference*- 1. 

When //i, M lf H 2 y and M 2 are set equal to zero, it 
will be seen that the resulting equations for stress 
are not the membrane stresses. This is true 
because tin* edge support is not the same as that 
necessary for the membrane condition when there 
are no edge shears or moments. This must be 
borne in mind if one wishes to separate the 
bending sti esses from the total stresses at any 
given point. 

The locations of the stress peaks were not 
determined for this case because of the complicated 
nature of the stress equations. 

Portion of a sphere loaded by edge shear and 
moment. — The following equations are applicable 
at any angle tp from the loaded edge of the spherical 
shell. The type and location of stress, as well as 
the rotation and deflection, are given by equations 
(42) to (45). Where ± or =F signs occur, the 
upper sign refers to the stresses on the inner 
surface and the lower sign to the stresses on the 
outer surface. See figure 12 for the sign conven- 
tion and figure 10 for curves of the functions 0, <£, 
and 12. Values of these functions are also 
tabulated in references 3 and 4. 

Meridional stress: 

|j^ cot *(*)]% 


Circumferential stress: 

*(M *(M] M, 


0(Xp)±-jj^y S2(X^)J // (l sin a (43) 


Rotation of meridian: 


2\ 2 T 2X ”1 

T — ^ O(Xv?)M 0 +^>(X^)/7 0 sin a\ (44) 


Outward displacement: 

2X 

Eh 


_ 2X 

w== Fl sin (a— 7?0(X^)// O sin a] 


(45) 

At the loaded edge, the previous equations 
reduce to 

6M 0 £Zo cos a 


(o*)*« o = dr- 


A* 


(40) 


/ ^ /2X 2 6A, r 2X tj . t . 

(<r*)*«o=( )M 0 — j-H 0 sin a (47) 

Vo= Eh(~R M 0 +H 0 sin a) =ut, 0 M 0 +u Ho Ho 


2X 

W<J= Eh sin a CM)— RTh sin a) = & Mo M 0 +& H(j H 0 


where 


(48) 

h 

(4 ) 


4X 3 2X 2 sin a 

Utr ° REh 


r _2X 2 sin a r 

Eh S "°— 


« Eh 

2 \R sin* a 


Eh 


(50) 


— jj- cot (a— ^)'P(Xp)±j^ 12(X<p)J H 0 sin a 


(42) 


The locations of the peak stresses are given by 
the following equations: 

Meridional stress: 



CSC 2 (a — <p m ) — 2X cot (a— <fi m )±6 j 

RH 0 sin a-t-2X 2 cot (a — </>„) M„ "j 

l^CSC 2 (a— <p m )±6j 

Rff 0 s'ma—2\ 

CSC 

* (a— <Pm)— x cot (a— <pj ± 6 J- 

mA 


(51) 


Circumferential stress: 

1 (3vjTX 2 )7?tf 0 sina±2X 3 M 0 

*Pq x ^ 7 I? v — - 

hp j±\ 2 )RH 0 sin «-6*X^M 0 


( 52 ) 


It is apparent that equation (51) is transcendental 
and must be solved by an iterative procedure. 
Because of this fact and the length of the equation, 
it might be more advantageous to plot the merid- 
ional stress distribution for a short distance using 
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equation (42). Then the location of the maximum 
stress, as well as its magnitude and the stress 
distribution curve, is obtained. 

The solution for the edge-loaded hemisphere 
can be obtained by simply substituting a=ir/2 in 
the previous equations. 

A discussion of edge influence coefficients for 
thin spherical shells is presented in reference 5. 
Three different approximate solutions, including 
the Geckeler approximation used in this section, 
are obtained for the differentia] equation of the 
thin spherical shell. The accuracy of these solu- 
tions and recommendations for their use are dis- 
cussed. 

Simply supported circular flat plate with uni- 
form load, edge moment, and edge normal load. — 

Refer to figure 7 for the geometry, loading, and 
sign convention for this case. The maximum and 
minimum stresses are found on the inner and outer 
surfaces at the center of the plate and are given by 


/ x , \ Ho, 0Mn T 3(3-Mq 2 

(<r r ) r *o= Wr-o* — jydr =F ^2 V (53) 


Sh 2 


where the upper and lower signs refer to the 
stresses acting on the inner and outer surfaces, 
respectively. The rotation of the edge of the 
plate is given by 

V== ~ (l+v)^ Mffl+ 8(l (54) 

The radial deflection of the edge of the plate is 

w=-£=^H 0 = 5„H 0 (55) 


Given or implied in equations (54) and (55) (see 
appendix D) is 


(1 +p)D 

uh= 0 


6m =0 


(1 —v)a 
Eh 


V (56) 


Z %il + v)D 


6 P = 0 


DISCUSSION 

In all shell junctions analyzed in this report, 
the shell is assumed to be long enough that one 
edge or junction can be analyzed independently 
of any other. In general, it is not necessary to 
make this assumption in order to solve for the 
induced forces at the junctions or for the stresses 
within a given shell. If the theory is available or 
can be derived for a given shell, it is only necessary 
to add two additional simultaneous equations for 
each junction to account for the close proximity of 
one junction to another. Then all equations for 
the junctions that are near each other are solved 
simultaneously. However, for most shells gen- 
erally used in engineering practice, the junctions 
are usually far enough apart to consider each one 
individually. Thus no serious limitation is in- 
curred by considering each junction separately as 
is done in this report. The work is considerably 
simplified because the mathematical difficulty in- 
creases very rapidly with the number of simultane- 
ous equations. 

The solution for the right circular cylinder is 
based on theory given in reference 3, as is that for 
the portion of a sphere. The analysis of the por- 
tion of a sphere is not valid in the region of the 
poles. For edge-loaded spherical segments with 
small included angles, shallow-shell theory such 
as in references 5 and 6 should be used. 

The analysis chosen for the cone (ref. 1) is 
essentially an “exact” analysis for relatively long 
cones, that is, those in which the edge effects do 
not overlap each other. In reference 7, by use of 
essentially the same basic differential equations, 
a solution is obtained for the short conical frustum 
with edge loads. The equations for both solutions 
are in terms of Bessel-Kelvin functions. Although 
it is not possible to present the results for the cone 
in as concise a form as some of the other shell 
solutions, most of the work has already been per- 
formed in putting them into tlieir present form. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, January 5, 1961 



APPENDIX A 

RIGHT CIRCULAR CYLINDER LOADED BY EDGE SHEAR AND MOMENT 


Consider a thin right circular cylindrical shell, 
subjected to axially symmetric shear forces and 
bending moments uniformly distributed along one 
edge, as in figure 9(a). Because of this loading, 
the internal forces and moments shown in figure 
9(b) arise. According to reference 3 (pp. 466 to 
471), the internal forces and moments are given by 

2V,= 0 

N 9 =2a&[m<p{px) +QoO(fa)] 

M x =Mj H0x) Q 0 Q(M 


Me = vM z — vM 0 $ Q 0 Q(fix) 

Q z =-2pM a n(fc)+Qo*(fa) 


The equations for the rotation V and the deflection 
w are given as equations (28) and (29) in the text. 
The functions 9, <£, and Q are shown graphically 
in figure 10. They are also tabulated in references 
3 and 4, The largest stresses at a given section 
corresponding to these forces and moments are 


ff ‘ v ’ = T [frV>*(/fr)+eoeGfc)l 
i7 'f 9 = j^A/ 0 '{>;/?/ ) Q o a(0.r)J 


Y (Al) 


The combined stresses in the axial and circumfer- 
ential directions are 


<r r = ± <Tm x (A2) 

and 

^=<^±<7^ (A3) 

respectively, and, with proper substitution, can be 
written as equations (26) and (27) in the text. 
In all instances where ± or T signs occur, the 
upper sign refers to the inner surface and the lower 
sign to the outer surface. 

Equations (26), (27), and (Al) give the stress 
distribution throughout the cylinder. In addition 
to the stress distribution, the locations and magni- 
tudes of the various maximum stresses are usually 
of considerable interest to the designer. The lo- 
cations where the various combined stresses are a 
maximum (except the possibility at the loaded 
edge) may be obtained by the usual theory of 
maximums and minimums (see the section 
METHOD OF ANALYSIS). Equating the de- 
rivative of r r with respect to x to zero gives 

x '4 an ' (A4) 

where x T is the distance from the loaded end of the 
cylinder to the location of the shear-stress peak. 
Similarly, the locations of the meridional and cir- 
cumferential stress peaks are obtained and given 
as equations (31) and (32) in the text. 


I!) 



APPENDIX B 


FRUSTUM OF A CONE LOADED BY EDGE SHEAR AND MOMENT AND INTERNAL PRESSURE 


The equations for the internal forces in a conical 
frustum subjected to edge shears and moments 
and internal pressure as shown in figure 11(a) can 
be obtained from references 1 and 8. These are 
given in terms of the Bessel -Kelvin functions as 
tabulated in reference 2, With the sign conven- 
tion as shown in figure ll'b), they ore 

Qu = ~~~ (abertS+abeiit 

+ C-i ker -2 £-C keA 2 £) (Bio) 


TL 


CSC a 

'~Y~ 


(Ci ber 2 Z+C 2 bei 2 Z+C 3 ker 2 % 


+ c ^ f) +fSh? (Blb) 


and the primes denote differentiation with respect 
to £. The equations for the rotation V and the 
deflection w are given as equations (35) and (36) 
in the text. The expressions for the constants 
of integration C\, C 2 , and C% are presented in 
the text as equations (37) for the loading condition 
shown in figure 11(a). 

The four constants of integration are deter- 
mined from the following boundary conditions: 

at £=£i=2mvVi 

H v =—H } My=M\ 

at f (B3) 

H v =-H 2 M v =M 2j 


A V = - (Ci ber 2 t+C 2 bei 2 £+C 3 ker 2 £ 

+ C t kei 2 & + Vy[ ^ a (Blc) 

Nt—Jj- (Ci ber^-^rCi bei'iC' C^kfi^i, 

& y 

+C A kei' t 1i)-\-py tana (Bid) 
[C,(£ bei'A+ 2v beiti) 

— C 2 (£ ber' 2 £+2v &er 2 £) + Cj(£ kei'^+ 2v kei 2 $) 

— C(£ ker'i £ + 2r ker 2 £) ] + (Bl(‘) 


Mt=— 2 7 ^ s \Ci( 2 bei 2 £+v£ beu£) 

— Ci(2 6cr 2 £+r£ bet’oO -\-Cz(2 kei 2 ^-\-y^ kei 2 £) 
-c t ( 2 iv/- 2 £+,£ ^£)]+ ^ (1 ^ } a (Blf) 


where 

£=2 M Vy 

m 4 ==12(l — p 2 ) ^ 
m 2 

** h tan a j 


(B2) 


As discussed in reference 3 (p. 563), the loads 
applied at one edge do not appreciably affect the 
stresses and deformations at the other edge if the 
cone is sufficiently long. It will be assumed that 
the thin shells considered here are long enough to 
conform to the above assumption. This is almost 
universally true for the conical frustums used in 
missile and spacecraft applications. 

The her and bei functions and their derivatives 
increase rapidly in an oscillatory manner as the 
distance y increases. Conversely, the her and 
kei functions and their derivatives which also have 
an oscillatory character decrease rapidly as the 
distance y increases. Hence, the terms involving 
her and kei and their derivatives are neglected 
when working near the large end of a conical 
frustum (or with a complete cone). The con- 
stants Ci and C 2 (which are associated with the 
her and bei functions and their derivatives) are 
then determined from the boundary conditions 
at the edge y=yi- Only when working with a 
frustum of a conical shell are all four constants 
of integration necessary. The constants C z and 
C 4 (associated with the ker and kei functions and 
their derivatives) are determined from the bound- 
ary conditions imposed upon the small opening of 
the truncated cone. 

For a long thin -walled conical shell whose half- 
angle a is not close to x/2, the calculations can 
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be simplified by solving two sets of two simulta- 
neous equations instead of one set of four equa- 
tions. C 3 and C 4 are neglected when evaluating 
Ci and C 2 from the conditions at the large end 
of the frustum, and C x and C 2 are neglected when 
evaluating C 3 and C 4 from the conditions at the 
small end. 

Substituting the boundary conditions at y=y x 
into equations (Bib) and (Ble) and neglecting the 
terms containing C 3 and C 4 , the following equa- 
tions are obtained : 


ber^+Ct bei&)- 

y i 


p ?/i sin 2 a 
2 cos a 


M i= 


2 m 2 yi 


[Ci(fi bei'Ji\-\-2v bei 2 i i) — U 2 (£i ber'<£\ 


+2v ber £ ,)]■ 


ph 2 tan 2 a 
8(1— v) 


Ci and C 2 are obtained from these equations and 
given as equations (37a) and (37b) in the text. 
Substituting the boundary conditions at y—y 2 
into equations (Bib) and (Ble) and neglecting the 
terms containing C x and C 2 yield the following two 
equations : 


-H 2 =— (C 3 ker 2 h+c 4 kei&) 
2/2 


py 2 sin 2 q 
2 cos a 


M 2 


-h 
2 m 2 y 2 


[C 3 ({ 2 kei' 2 Z 2+2 v keioQ — C^z ker^ 2 


+2v &er 2 £ 2 )H 


ph 2 tan 2 a 
8(1 — v) 


From these equations, C 3 and C 4 are found and 
presented as equations (37c) and (37d) in the text. 

The closed-form solutions for the four constants 
of integration can now be substituted into equa- 
tions (Bl) to determine the complete distribution 
of internal forces. As was true in the derivation 
of the constants, the terms containing C 3 and C 4 
can be neglected when working near the large 
end of a long frustum (or on a complete cone), 


and the terms containing Ci and C 2 can be neg- 
lected when working near the small end of a 
frustum of a long thin-walled conical shell. 

The influence coefficients for the edge rotation 
and deflection of a truncated cone can now be 
found. Substitute the expressions for the con- 
stants of integration (eqs. (37)) into equations 
(35) and (36) and solve for the edge rotations and 
deflections. The coefficients of the applied loads 
wdll be the desired edge influence coefficients, as 
shown symbolically in equations (35a), (35b), 
(36a), and (36b). 

The relations between edge loading and edge 
rotation and deflection can be expressed in terms 
of the edge coefficients as equations (35a), (35b), 
(36a), and (36b) in the text, where the subscript 1 
refers to the edge of the large end (y—y i) and the 
subscript 2 to the edge of the small end (?/ = 2 / 2 ) . 
Using the loading condition shown in figure 11(a) 
and the sign convention of figure 11(b), the expres- 
sions are obtained for the edge influence coeffi- 
cients and given as equations (38) in the text. 

The meridional stresses and the circumferential 
stresses are given by 


Ny 6 My 
h ± h> 

(B4) 

JV$ ( 6il /, e 

(B5) 

h ± v 


respectively. With the proper substitutions, these 
equations become equations (33) and (34) in the 
text. In all instances where ± or T signs occur, 
the upper sign refers to the inner surface and the 
lower sign to the outer surface. At any given 
section the maximum shear stress acting on a 
middle surface is given by 


3 Q y 3 cot a 
2 2 hy 


(Ci ber^+C-z bei 2 £ 


+C 3 kerg+Ct keiit) (B6) 


Equations (33), (34), and (B6) give the complete 
stress distribution throughout the conical frustum. 



APPENDIX C 

PORTION OF A SPHERE LOADED BY EDGE SHEAR AND MOMENT 


A portion of a sphere subjected to rotationally 
symmetric uniformly distributed edge loading, as 
in figure 12(a), will have present internal forces 
and moments shown in figure 12(b). According 
to reference 3 (pp. 549 to 551), the internal forces 
and moments are given by the following (with a 
change in notation for <p ) : 

Av=— J^3/oQ(X#>)+flo (sin a) ^(X^Jcot («—*>) 

9\2 

A 7 * M o* (\<p) — 2 \H 0 (si n a) 0 (\<p) 

H 0 (sin a)Q(\<f>) 

A 

Mt=vM t =vMo$(\<p)— ~ H 0 (sin a) Q(\<p) 

A 

o\ 

M 0 Sl(\<p) +Ho (sin 


The equations for the rotation I 7 and the deflection 
w are given as equations (44) and (45) in the text. 
The functions 0, 4>, and ft may be obtained 

from figure 10 or the tables in references 3 and 4. 

In this analysis, circumferential stresses (sub- 
script 0) refer to the stresses normal to the cross 
section obtained by cutting the spherical segment 
with a plane containing the axis of revolution. 
Meridional stresses (subscript <p) are those acting 
on the cross section cut by a cone whose apex is at 
the center of the sphere and whose axis coincides 
with the axis of revolution of the spherical seg- 
ment. 

The largest stresses at a given section corre- 
sponding to the various internal forces and mo- 
ments are 


<rv,=^=- COt( ;* y) [|M,O(M+g 0 (sin «)*(X*o]' 

Ho (sin a) e(M 
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(Cl) 


The combined stresses in the meridional and 
circumferential directions are 


and 




(C2) 

(C3) 


respectively, and can be written (with proper 
substitution) as equations (42) and (43) in the 
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text. Where ± or T signs occur, the upper sign 
refers to the inner surface and the lower sign to 
the outer surface. 

Equations (42), (43), and (Cl) present a com- 
plete picture of the stress distribution throughout 
the spherical segment. From these equations, 
the locations and magnitudes of the various 
maximum stresses can be obtained as a further 
aid to the designer. The points wdiere the com- 
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bined stresses are a maximum (except for the 
possibility at the loaded edge) may be found by 
the usual theory of maximums and minimums 
(see the section METHOD OF ANALYSIS). 
Setting the derivative of r+ with respect to <p 
equal to zero and solving for <p give 

1 /. R Hq sin oA 

arctan^l-y ^ J (C4) 

where <p T is the angle between a radius at the 
loaded edge and a radius (in the same meridional 


plane) at the location of the shear-stress peak. 
The locations of the meridional and circumfer- 
ential stress peaks, given as equations (51) and 
(52) in the text, are obtained in a similar manner. 

The equations presented in this appendix apply 
to a portion of a spherical shell whose half-angle 
ot may be less than or greater than 7r/2 but does 
not approach 0 or ir (fig. 12). The reason for 
this is that shallow-shell theory (refs. 5 and 6) 
must be used near the poles. The special case of 
the edge-loaded hemisphere may be obtained by 
substituting a=x/2. 


APPENDIX D 


SIMPLY SUPPORTED CIRCULAR FLAT PLATE WITH UNIFORM LOAD, EDGE MOMENT, AND EDGE 

NORMAL LOAD 


Consider 1 lie ease of a simply supported circular 
flat plate under the action of a uniformly distrib- 
uted normal load p and with a radially symmetric 
edge normal load H Q and an edge moment M 0 as 
shown in figure 7. The plate can most readily be 
analyzed as three separate problems and then, by 
superposition, the final results can be obtained. 

Consider first the action of H 0 on the plate. 
This is a case of hydrostatic plane stress: 



From reference 3 (p. 43) it can be shown that the 
rotation of the edge of the plate due to M 0 is 


V 


aM 0 
(1 -\-v)D 


=&m q M o 


Now consider the simply supported plate acted 
upon by the uniformly distributed load p. From 
reference 3 (p. 57) the maximum stresses occur at 
the center of the plate and are shown to be 


(y^r)mat — (y fl) max — “F 


3(3+v)a* 
Sh 2 


The radial deflection at the edge is 

The surface of the plate remains flat, and therefore 

Wtf 0 =0 


For small deflections normal to the plate, the 
radial deflection can once again be assumed zero. 
Therefore, 

S p =0 

The angle of rotation of the edge can be shown 
to be 


Next analyze the plate under the influence of 
M 0 . The moment at any point in the plate is M 0 . 
Therefore, the stresses at the inner and outer sur- 
faces are 


For small deflections normal to the plate, the radial 
deflection is zero and 

K=o 


V ~~8(l+v)D p ~ WpP 

By adding the effects of the conditions described 
previously, the following results are obtained . The 
maximum and minimum stresses are found on the 
inner and outer surfaces at the center of the plate 
and are given by equation (53) in the text. The 
edge rotation and deflection are given by equa- 
tions (54) and (55) in the text, respectively. 







APPENDIX E 

SHEAR AND MOMENT AT AN AXIAL CHANGE OF THICKNESS IN A CIRCULAR CYLINDER 


In this appendix and in the following appendixes 
the discontinuity shear force and bending moment 
are determined for several specific shell junctions. 
Because the derivations for the shear force and 
bending moment at the various junctions are very 
similar, they are given here in a general way. 
These equations are then used with the proper 
subscripts for the shear and moment in most of 
the following junction problems. 

The general equations for the so-called edge 
rotations and deflections of two shells, i and j , at 
their junction can be written in terms of internal 
pressure, discontinuity shear force and moment, 
and edge influence coefficients as 


w 0t t — oH - pp 

Vo, i=o) iiHo Ho+(ai.M 0 MQ+(a it P p 

(h.1; 

5j.h q Ho+5j.m 0 M 0 + 8j, P p 

Vo t j = WjM Q HQ+Wj,M Q Mo + Wj, p p J 

The equations for continuity of rotation and de- 
flection at the junction are 

m 

Wo,i= w o,j J 

By substituting equations (El) into equations 
(E2), the unknown discontinuity shear force and 
moment are found to be 


° (5< f jr 0 “ &) tf 0 ) («f.Af 0 “ w^a/ 0 ) — «y.ff 0 ) 


These equations can now be used with the proper 
subscripts and edge coefficients to determine the 
discontinuity forces at most of the junctions 
considered in this report. 

The case considered in this appendix is that of 
a radially symmetric longitudinal change in the 
wall thickness of a thin-walled circular cylindrical 
pressure vessel. The radial expansion of a pres- 
surized cylinder with closed ends due to the mem- 
brane stresses is 

ms ' p iE5) 

This equation applies only at a distance from the 
closed (Uids where the effects of the closures are 
negligible. The membrane stresses cause no ro- 


tation of the meridians of a cylinder, and therefore 

u>p=0 (E6) 

The free radial expansion of a cylinder due to inter- 
nal pressure is inversely proportional to the wall 
thickness. Consequently, discontinuity forces 
and moments must be present which will make the 
displacements and rotations of two adjacent re- 
gions of different wall thicknesses identical at 
their junction. 

Assume the discontinuity forces to be directed 
as shown in figure 1(a). Refer to equations (E3) 
and (E4); for this case, H n = Qo- Letting sub- 
scripts 1 and 2 refer to the regions of thicknesses 
ht and h 2 , respectively, the expressions for the un- 
known discontinuity shear and moment are (with 
w p =0) 


— (5i, g — 1 (E7) 

* " _ (5 1 , <J 0 — 5a , o 0 ) (oi 1 , — w 2 . M 0 ) — (Si . -Vf 0 — S 2 , m q ) Oi . <3 0 — "2. 0 0 ) 

A[ (Si.p— frj.pXm.QQ— M2.o 0 ) _ P ( Es ) 

(5 1 . <? 0 — 5 2 , « 0 ) ( W i , ,w 0 — M 2, .V„) — (5l , M 0 — ,M a ) (« I . O 0 _ 
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11 the discontinuity shear and moment are di- the influence coefficients obtained from equations 

reeled as shown and the pressure is internal, then (E5), (E6), and (30) are as follows: 


* (l v \ q2 
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^2 .*r 0 = 
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Substituting the edge influence coefficients as chosen in figure 1(a), the sign for Q 0 must be 

given in equations (E9) into equations (E7) and changed when solving for the stresses in the cyl- 

(E8) results in equations (1) and (2) in the text for inder of thickness h 2 . 

Qo and M 0 . Because of the sign convention 



APPENDIX F 

SHEAR AND MOMENT AT AN AXIAL CHANGE OF THICKNESS IN A CONE 


The radial expansion of a conical shell due to 
internal pressure is a function of the thickness; 
therefore, if the thickness varies along a meridian, 
discontinuity forces are induced to make the slopes 
and deflections continuous along the meridian. 
The discontinuity shear and moment are deter- 
mined here in terms of the edge influence coeffi- 
cients and internal pressure for a truncated 
conical shell with a radially symmetric sudden 
change of thickness at a distance y 0 from the apex. 
Positive shear and moment are as shown in figure 
2 . 

Letting the subscripts 1 and 2 refer to the regions 
of thicknesses h x and h 2f respectively, the expres- 
sions for the discontinuity shear and moment can 
be written as equations (4) and (5) in the text by 
referring to equations (E3) and (E4). 

The expressions for the edge influence coeffi- 
cients for the cone are given in the section “Frus- 


tum of a cone loaded by edge shear and moment 
and internal pressure.” For the problem being 
considered here, 2 /i=2/2=2/o and r 1 = r 2 =r 0 , but X 
and £ are not the same for the two edges at the 
junction because hi does not equal h 2 . The 
parameters \i, £i, «i, and Si are associated with the 
region of thickness hi ; likewise, X 2 , £ 2 , and 5 2 
are associated with h 2 . The expressions for the 
edge influence coefficients necessary for the solu- 
tion of H 0 and M 0 are given explicitly as equations 
(7) in the text. 

The desired discontinuity shear and moment 
are found by use of equations (4) and (5), respec- 
tively. The deflection, rotation, and internal 
stresses at any point in the cone can now be de- 
termined by use of the equations in the section 
“Frustum of a cone loaded by edge shear and 
moment and internal pressure.” 
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APPENDIX G 

SHEAR AND MOMENT AT A CHANGE OF THICKNESS IN A PORTION OF A SPHERE 


Consider the case where a change of thickness 
takes place in a thin-walled spherical shell, sym- 
metrical about an axis as shown in figure 4. The 
radial expansion due to the membrane stresses 
only is different in the regions of different thick- 
nesses. The shear and moment required to make 
the slopes and deflections of the two regions 


coincide at their junction are determined in the 
following paragraphs. 

Let the subscripts 1 and 2 refer to the regions 
of thicknesses h\ and h 2 , respectively. Then the 
expressions for the discontinuity shear and moment 
can be writ ten in the following manner by referring 
to equations (E3) and (E4) : 


1 , H 0 &2, H 0 ) («1 , M 0 &2, A/ 0 ) — 1 , M 0 — ^2, M Q ) (&>1 , — Ci> 2 , Q ) 


M 0 = 


(5].p — $2 — W 2 ,H() 


— <^2.AT 0 ) — (*5l,Jtf 0 5 2 ,A/ 0 )(c*>i i ff 0 — £02,^) 


V 


(Gl) 

(G2) 


where £*>i,p=w 2 ,p=0. It is known from con- 
ventional membrane theory that the displace- 
ment perpendicular to the axis of symmetry of a 
point on the shell due to membrane forces only is 

2 £7jT OL = ^vP (G3a) 

Since the membrane forces cause no rotation 
of the meridians, 

cd p =0 (G3b) 

The other necessary edge influence coefficients 
can be found in equation (50). 

With internal pressure and the discontinuity 
shear and moment in the directions shown, the 
edge influence coefficients are as given in the 
following equations: 
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_ 4X 3 2 

Ul,A/ 0 = 

REh i 

W 2 ,.V 0 = 

~REh 2 J 


If the influence coefficients of equations (G4) 
are substituted into equations (Gl) and (G2), the 
equations for 77 0 and M 0 take the form of equa- 
tions (11) and (12) in the text. 


APPENDIX H 

SHEAR AND MOMENT AT THE JUNCTION OF A CYLINDER AND A CONE 


In missile or space structures, a conical dome 
or bulkhead is usually attached to the main 
cylindrical body by means of a transition section 
such as a torus. However, there are applications 
where the cone is welded directly to the cylinder 
(fig. 5). For such cases, the discontinuity shear 
and moment at the junction are determined in 
the following paragraphs. 

Let the subscript c refer to the cylinder and k 
to the cone. Positive shear and moment are as 
indicated in figure 5. From equations (E3) 
and (E4), the discontinuity shear and momeni 
are determined as equations (14) and (15) in 
the text. 

From membrane theory, the edge-displacement 
influence coefficient for the cylinder is 



Since the meridian of a cylinder does not rotate 
because of membrane stresses, die edge-rotation 
influence coefficient is 

w C(P = 0 (H2) 

The other necessary edge influence coefficients 
are given in equations (30) and (38), where the 
coefficients for the large end of the cone (subscript 

1) are used. There are instances when the 
cylinder may be attached to the small end of the 
conical frustum. For this case, the coefficients 
for the small end of the conical frustum (subscript 

2) should be used. 
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APPENDIX I 


SHEAR AND MOMENT AT THE JUNCTION OF A CYLINDER AND A PORTION OF A SPHERE 


The cylindrical pressure vessel is frequently 
closed at the ends by a portion of a sphere. The 
dome may subtend an angle considerably less 
than 180° (fig. 6(b)). The special case of a 
hemispherical dome as the closure (fig. 6(c)) is 
commonly used. The equations given here are 
also applicable to the case of a subtended angle 
in the dome greater than ISO 0 (fig. 6(d)). The 
restrictions imposed upon the spherical shell in 
appendix C appl}~ here also. The discontinuity 
shear and moment at the junction of the cylinder 
and spherical dome are determined in the following 
paragraphs. Positive shear and moment are as 
shown in figure 6(a). 

With the subscript c referring to the cylinder 
and s to the portion of a sphere, the following 
equations express the edge rotation and deflection 


of the two si jells at their junction : 


w 0 ,r 5 cMq H 0 +5 cMo M 0 +5 c . pP 

V o, c=<*>c,n Q Ho J ru c ,M 0 M 0 J ru) Ct P p 



Continuity of rotation and deflection at the 
junction requires 


V 0 , e =V 0t9 

tt'o, 0, * 


} 


( 12 ) 


Substituting equations (II) into (12) gives 


H, 


0— ^ ,M 0 ) — ^5.// 0 ( W C,A/ 0 ““W5.3/ 0 )] 2 COS a 

+ — ($ Ct A f{) — $,,A/ 0 ) — Vs.p) ^ V 

"*"[(*■ C,A/ 0 — 5 s ,M 0 ) U*,H q ) — (Sc,H 0 — &s,H q ) (fc>c,Af 0 — Us,M q )] 
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7? 


(Ss.HqUc'Hq 9~ COS C¥ + (&e,H Q °Vp) (&c. p &t, p) fac,H Q 

£ — _ . , ** — p 

(<5 c,M 0 5, > ,\f 0 )(aJ C(iWo W Si ff 0 ) (dc,H 0 &»,H 0 )(Vc,M 0 Vs.Hq) 


(13) 

(14) 


The edge coefficients for the membrane stresses 
are 

/ A a* (1 — v)R* sin a 

2.fe 5s ' p ~ 2EK 

w fiP =0 =0 
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The remaining edge influence coefficients are 
obtained from equations (30) and (50). Substitu- 
tion of the edge coefficients into equations (13) 
and (14) results in equations (17) and (18) in 
the text. 




APPENDIX J 

SHEAR AND MOMENT AT THE JUNCTION OF A CYLINDER AND A FLAT HEAD 


One of the simplest closures for a cylindrical 
pressure vessel is the flat plate head. This is a 
very stiff restraint on the cylinder and therefore 
induces stresses in the cylinder higher than would 
be produced by most of the other closures. 
However, since it is the easiest and cheapest to 
fabricate, it is sometimes used in experimental 
facilities or for other such purposes where the 
weight of the vessel is not critical or where the 
higher induced stresses can be tolerated. 

The discontinuity shear and moment are deter- 
mined here for this case with the sign convention 
being as shown in figure 7. Let the subscript c 
refer to the cylinder and / to the flat plate. With 

w c ,p=5 / ,3f 0 =5 /t p=w /iHo — 0 (Jl) 


equations (E3) and (E4) become 

rr — p — be.JUc. M (i W/. M( ) „ , T o\ 

H 0 =-7Z : r? s J P 

( Oc , H Q — Of. H 0 ) (“c. Jlf Q — 01 f. M 0 ) ~0 e . i f Q U c . H 0 

and 

,, K pWc. H 0 +u r . , (S c , B — 5 f , H 0 ) 

M 0 =t7 : — w s — i v (Jo; 

(oc.H 0 -~Of t H Q )(U Ci M 0 Oc,M q Wc,Hq 


Substituting the expressions for the edge 
influence coefficients found in equations (16), 
(30), and (56) into equations (J2) and (J3) yields 
equations (20) and (21) in the text. Having found 
the desired discontinuity shear and moment, the 
stresses in the cylinder and plate can now be 
computed. 


31 



APPENDIX K 

SHEAR AND MOMENT AT THE JUNCTION OF A CONE AND A PORTION OF A SPHERE 


The toriconical head is frequently used in tank 
design. The special case in which the torus is a 
portion of a sphere (fig. 33) can be analyzed in 
the following manner. Because the expansions of 
the cone and the portion of a sphere due to mem- 
brane stresses only are not, in general, the same 
at the junction, discontinuity forces must be 
present to make the rotations and deflections 
identical at their common boundary. 

Assume positive shear and moment to be in the 
directions shown in figure S. The rotation and 
deflection at the junction of the cone and the 
portion of a sphere can be expressed in terms of 
edge influence coefficients, discontinuity shear and 
moment, and internal pressure. The equations 
expressing the relations can be written as follows: 

W 0 , * = $k, * 0 Ho *T &k, MqMq -j- S k ' pp 

MqMq + pp 

= COS + vP 

ko, s = U$,H 0 cos ^ ~f* <*>s, pP 

(Kl) 

The subscripts k and s refer to the cone and spheri- 
cal shell, respectively. Continuity of rotation 
and deflection at the junction results in the 
following equations: 
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Wo f fc=W 0> , 


} 


|(K2) 


Substitution of equations (Kl) into equations 
(K2) gives equations (23) and (24) in the text. 
The edge influence coefficients are obtained from 
equations (38) and (50). Note the changes in 
the signs of some influence coefficients, as men- 
tioned in the body of the report. In addition, the 
coefficient is zero, since the meridians of a 
sphere do not rotate when subjected to membrane 
stresses only. The edge-displacement coefficient 
for the portion of the sphere due to internal pres- 
sure d 8tP can readily be obtained from membrane 
theory. The expressions for co StP and & S , P are 
then 



Figure 13. — Toriconical head on cylinder. 





APPENDIX L 

EFFECT OF NONCONCURRENCE OF MIDDLE SURFACES AT SHELL JUNCTIONS 


Frequently the middle surfaces of two shells do 
not coincide at their junction. This situation 
may arise because the shell has been milled on 
only one surface or because two shells of different 
thicknesses have been joined to give a smooth 
outer contour or for some other reason. 
This eccentricity of the middle surfaces induces 
bending stresses in addition to those previously 
found for the shell junctions which have continuous 
middle surfaces. The following paragraphs show 


H 



Axis of 


revolution-^ 
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Z 


Axis of revolution^ 

(b) 

(a) Fundamental Force system. 

(b) Conventional force system. 

Figure 14. — Shell junction with nonconcurrent middle 
surfaces. 


how these additional bending stresses can be 
found. 

Consider first the case where the tangents to 
the meridians at the shell junction are parallel. 
Let the meridional membrane force be denoted by 
N and the distance normal to the shell between 
the middle surfaces be d as shown in figure 14(a). 
As before, imagine the shell to be cut into two 
pieces by a plane passing through the junction. 
Very fundamentally, it is possible to represent the 
forces at the cut by equal and opposite shear 
forces H and equal, opposite, and colinear normal 
forces N acting along some unknown line of action. 
Assume this line of action is at distances / and g 
from the middle surfaces of shells i and j, respec- 
tively. Since N represents the resultant meridi- 
onal force on either side of the cut, there is no 
bending moment present. The meridional forces 
are moved to the middle surfaces along with 
bending moments NJ and Ng as shown. Notice 
that the moments acting on the two edges are not 
equal. This is a result of the nonconcurrence of 
the middle surfaces and the fact that the meri- 
dional forces N are not now colinear. Let 

M t =Nf (LI) 


Then, considering only the geometry, the moment 
on the other side of the junction is 

M j =Ng=M-Nd (L2) 


The equations of continuity at the junction are 

Wo,i=WQ'jJ 

In terms of edge influence coefficients, edge loads 
and internal pressure, equation (Lo) becomes 

WlffH i + Wf, pP=Uj t + 0)j j -T p p 

j+^j.pP 

(L4) 
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Substituting equation (L2) into (L4) and solving for JT and M t yield 

JJ W yj|/) ($i, p ^ j, p) (w^ p Wy 7 ?) *5 j ,Af) ^ 

(^i,H W; ,h) {&i,M & j.Af) &j.Af) (&i,H &j.h) 


W ; . Af 5 j, Af 0?/, jf 5 j i a/ 


(w J( // Wj,//) (Si'M — Sj'M) — (Wi'M — £d>, m) (&i,H — bj.n) 


Nd (L3) 


/ fe/. u W;, bj.tt) (w,,g W>. //)(&<■ p p 

{&i,H — ^j.Af) — (<*>* , M — Wy, Af) ($ i . J7 — &j.h) 




(w j , H Wj.ff) (&i,Af &J,m) — Wy t m) ( 5 i , // — 5/, //) 


Atf (L6) 


Substitution of equation (L6) into (L2) gives 


^ __ (w,. p — w ; -. (wj,y — Wj.//) (hj' V — 6j, p ) 

J (w,.// Wj //) (6*,*/ — 5 Jfi v) — (Wf.il/ — itij' m) (5i'H — $j.h) 


V 


w,.ii/(5j H 8j,h) (wj l g^W/ L g)£i. 

— /.I. -A ^ — A. »A — (f.\. . ^ M , . A ^ } 


M 


A T /7 .71 


The first parts of equations (L5) to (L7) are 
recognized to be the discontinuity forces deter- 
mined before when the middle surfaces were con- 
tinuous (see eqs. (E3) and (E4)). The other 
terms are the effect of the nonconcurrence of the 
middle surfaces. Thus the discontunity forces 
due to only the nonconcurrence of the middle sur- 
faces can be computed by using the second half of 
equations (L5) to (L7). These forces can then be 
added to those given in the body of the report to 
obtain the resultant stress distribution. 

The sign of the bending moment Nd must be 
assigned with care. It would be oppositely di- 
rected if the middle surface of shell i were inside 
the middle surface of shell j at the junction (fig, 
14(b)). 

If the tangents to the meridians at the shell 
junction are not parallel, it is not possible to use 
the foregoing equations because the meridional 
forces are not equal and opposite. In such a case, 
it is probabl3’ most convenient to work with forces 
at the junction which are either parallel or per- 
pendicular to the shell axis. As in the previous 
analysis, the components perpendicular to the 
shell axis do not affect the edge moment. In com- 
puting the bending moments due to nonconcur- 
rence of the middle surfaces in this case, the eccen- 


tricity of the middle surfaces must be measured 
normal to the shell axis. Equations similar to 
(L5) to (L7) will finally result. Because of the 
lack of availability of all necessary influence coef- 
ficients, this case was not pursued further. 
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